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Abstract 

We study the properties of the nuclear rotational excitations with hypothetical tetrahedral 
symmetry by employing the microscopic mean-field and residual-interaction Hamiltonians with 
angular-momentum and parity projection method; we focus on the deformed nuclei with tetra- 
hedral doubly-closed shell configurations. We find that for pure tetrahedral deformation the ob- 
tained excitation patterns satisfy the characteristic features predicted by group-representation 
theory applied to the tetrahedral symmetry group. We find that a gradual transition from 
the approximately linear to the characteristic rigid-rotor, parabolic energy-vs.-spin dependence 
occurs as a function of the tetrahedral deformation parameter. The form of this transition is 
compared with the similar well-known transition in the case of quadrupole deformation. 

PACS numbers: 21.10.Re, 21.60.Jz, 23.20.Lv, 27.70.+q 
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I. INTRODUCTION 



Symmetries play an important role in physics, often serving as guide-lines in studying 
characteristic features of motion of quantum systems. In particular, in nuclear physics, the 
spatial symmetries of the nuclear mean-field potential are crucial in determining both the 
properties of the independent particle motion and of the nuclear stability. The best known, 
the spherical symmetry of the potential, implies the high degeneracy of the single-particle 
level energies, the so-called magnetic degeneracy, (2j + 1), of the orbitals characterised by 
the spherical-shell angular- momentum quantum- number, j, and leads to the well-known 
spherical magic numbers in nuclei. The spontaneous breaking of the spherical symmetry 
arises when the energies of the spherical configuration of the system with certain particle 
numbers are higher as compared to the energies of the alternative non-spherical spatial 
distribution of nucleons, the mechanism to which the Pauli exclusion principle contributes 
importantly. The spherical symmetry breaking removes the (2j + l)-degeneracy and leads 
to the deformed single-particle orbital scheme 

More generally, nuclei which are not spherical may acquire the forms governed by 
the point-group symmetries, some being more likely than the others - depending on 
the actual number of nucleons. In this context, it has been suggested, cf. Ref. and 
references therein, that each symmetry group 'sufficiently rich' in terms of symmetry 
elements may lead to its proper scheme of magic numbers, which in turn can be seen 
as characterising such a symmetry group from the point of view of realistic realisations 
of the nuclear mean-field theories. In particular, a series of earlier publications related 
to point-group symmetries focused on the tetrahedral and octahedral symmetry groups 
which are the only ones that lead to an extra (four-fold) degeneracy of single-nucleonic 
levels in deformed nuclei and to an increased nuclear stability. In other words: The 
tetrahedral and octahedral symmetries are the only ones within which both the four- 
fold degeneracy ('unusual' case) appears for some levels and the two- fold degeneracy 
('usual' case) appears for some other levels. This should be compared to the habitual 
two-fold (Kramers) degeneracy associated with all the levels known in all other deformed 
nuclei as e.g. in the case of the quadrupole deformation, see e.g., Ref. In fact, the 
symmetry-implied large shell gaps occur at some specific nucleon numbers and they can 
be comparable to the spherical shell-gaps. The tetrahedral magic numbers, Nt and Zt, 
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for the neutrons and protons, respectively, are {4]: Nt 01 Zt = 16, 20, 32, 40, 56, 68-70, 
90-94, 112, and 136/142. 

The possibihty that tetrahedral symmetry is present in atomic nuclei has been dis- 
cussed as early as in 1970'ies for ^^O in relation to the hypothetical four alpha-cluster 

91], or the 



structure 



5|-l8l]. Calculations employing the 



microscopic-macroscopic method 



Skyrme Hartree-Fock (HFB) method 10Hl2[ suggested that, in heavier nuclei, the tetra- 
hedral shapes may appear in low-l ying excited-, or even in the ground-states for specific 

3-161 refore, both timely 



nucleon numbers; see e.g. Refs. 



and interesting to employ the well established methods of the theory of nuclear structure 
in an attempt of examining the leading features of the excitation spectra of collective 
motion associated with of the tetrahedral shape. 

In this article we focus on certain properties of nuclear rotational bands in tetrahedral- 
deformed nuclei using advanced microscopic techniques which employ the angular- 
momentum and particle-number projection methods. Although the methods of perform- 
ing the projections are straightforward and well-known j^, 17], their numerical realisation 
is a non-trivial task, especially for the non-axially symmetric and non-time-reversal in- 
variant systems described using the cranking approximation (see e.g. Refs. 18|, ll9|). We 
have developed an efficient method to perform the projection from general and realistic 



mean- field wave functions calculated with large number of basis states [20|. Characteristic 
features of our method which have an important impact on an increased efficiency in the 
numerical realisation of the algorithm can be summarised as follows: (a) An efficient trun- 
cation scheme using the information about the occupation probabilities in the canonical 
basis; (b) The full use of the Thouless amplitude with respect to a Slater determinant 
state in place of the generalised Bogoliubov amplitudes; (c) Avoiding the sign problem 
for the norm overlap in terms of the Pfaffians [2l|] with using the Thouless amplitude. 

The tetrahedral-symmetry nuclear-states have not been so far identified in nature. In 
order to facilitate the associated experimental research program in the case of a possible 
discovery of a new quantum phenomenon, one needs to establish first of all the global and 
leading characteristic properties of the tetrahedrally-symmetric nuclei. The main focus 
of the present article is to test the projection techniques associated with the mean-field 
methods, as far as the characteristic features of the energy spectra are concerned rather 
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than trying to be as realistic as possible in terms of the energy scale predictions. Another 
aspect is to test the projection techniques necessary to calculate the electromagnetic 
transition probabilities within the nuclear-mean field theory. Such transition probabilities 
and/or their branching ratios can be used as characteristic signs of the symmetries and 
will become a necessary tool for establishing such symmetries in nature. The related 
research program is in progress and results will be published elsewhere. 

In the next Section, Sec. [Til the principal mathematical expressions of the method are 



briefiy recapitulated; the details can be found in Ref. [20|. The results of calculation are 
presented and discussed in Sec. Illli The final section, Sec. HVl is devoted to the summary 
and possible future perspectives. 

II. METHOD OF CALCULATION 

In what follows we assume that a nucleus is in a state corresponding to a tetrahedral- 
symmetry minimum, one possibly among other competing minima in the total energy 
landscape. Assuming purely static configurations, i.e., ignoring the collective effects such 
as the zero-point vibrations or any other form of, e.g., large-amplitude motion which may 
be particularly needed in the case of the flat energy landscapes, we will calculate the 
excitation pattern using the angular-momentum, parity, and particle-number projection 
techniques, the latter in relation to pairing, combined with the mean-fleld techniques. 

Thus, the most general symmetry-conserving wave function is sought in the form 

K 

where Pijx and P± are the angular-momentum and parity projectors (see e.g. Ref. 
The mean-fleld state 1$) is taken in the form of the anti-symmetrised product (HFB 
type) wave-functions, which is specified in Sec. Ill A[ The i^-mixing coefficients, g^^"^-, are 
obtained by solving the generalised eigen- value problem of the Hill- Wheeler equations 

K' K' 

with the Hamiltonian and norm kernel matrices being defined as usual as: 

I = (^1 I pkK'P±m- (3) 




In the present approach we wish to go beyond the mean-field approximation without 
perturbing the tetrahedral symmetry of the problem. This can be done by introducing in 
Eq. ([3]) a two-body spherically-symmetric Hamiltonian whose form will be discussed 
in Sec. HTBI 

The neutron and proton number projections require the number projectors (P^ and 
P^), which are further included in Eqs. (II]) — ([3]). However, we found that the effect of the 
number projections on the quantum spectra in the present work is small, see Sec. IIII Bt 
and they are simply neglected in most cases after verifying that such a neglect is justified. 



A. Mean-Field Model and Its Hamiltonian 



It is often of interest to employ the consistency condition between the mean-field and 
the many-body Hamiltonians like, e.g., in the Skyrme-Hartree-Fock models, since it is 
believed that such a consistency offers a more realistic description of the many-body 
systems in question. In this article we wish to focus first of all on the nuclear excitation 
spectra in the presence of pure tetrahedral symmetry of the mean-field. In this context it 
is preferable to work with the model allowing to completely control the deformation and 
the underlying geometrical symmetry. For this purpose, a phenomenological deformed 
mean-field is more convenient. 

We use the product-type states composed of the eigen-f unctions of the Woods-Saxon 



(WS) potential 22|, for which the general deformed shape is parametrised with the help 



of the spherical-harmonic {YA^t}-basis: 



(4) 



where the coefficient c^,({a}) takes care of the volume-conservation condition. In the 
present application, the pure tetrahedral deformation is realised by requiring all deforma- 
tion parameters = except 032 (more precisely, 03+2 and 03-2) • In this particular 
case, the problem of the centre of mass does not arise since tetrahedral-symmetric uniform- 
distributions preserve the position of the centre of mass independently of the value of the 
tetrahedral deformation. The deformed WS single-particle Hamiltonian, hdd, is diagonal- 
ized in the spherical harmonic oscillator basis with the oscillator quantum numbers n and 
/ satisfying the usual relations + Uy + Hz = 2n + I < 



' max • 



The HFB type product state is obtained by further including the monopole pairing 
field with the particle number constraint: 



h 



A^^PI + P^)- KNr. (5) 



pair / J 

T=n, p T=n, p 



The pairing gap A in the Hamiltonian above is either parametrised, cf. Sec. Ill C\ or self- 
consistently calculated by using the HFB treatment assuming the seniority interaction for 
neutrons and protons, 

Hp = -J2GrP!Pr- (6) 

T=n, p 

The simple-minded usage of this pairing interaction with large model space may lead to 
a divergence. We replace the monopole pairing operator, P\ by the one with a cut-off 
function, /c(e), 

pt = ^44 ^ ^/,(e,)44 (7) 

i>0 j>0 

where the quantity and d] are the eigen-energy of the deformed WS single-particle 
state and its creation operator, respectively, so that h^ci = J2i ^i^ldi, and i refers to the 



time-reversal conjugate-state of i. The form of the cut-off function is chosen to be 23|, 



m = \ 



1 + erf 



e — A + Alow 



d. 



cut 



1/2 



1 + erf 



-e + A + Kp 



d. 



cut 



1/2 



2 r _2 

with the error function defined as usual by erf(x) = I e ^ dt. Following Ref. 23 1 

V Jo 

the parameter values adopted here are: Aup = Aiow = 1.2 hw and dcut = 0.2 hw with 
huj = Al/A^/^ MeV. The chemical potential A in Eq. ([8]) is simply chosen as (ejg -|-ej(,+i)/2, 
where io is the last occupied orbital in the case of no pairing. Such a treatment results in 
preserving typically -|-(25-to-35) and — (15-to-25) states around the Fermi level for Rare 
Earth nuclei depending somewhat on the deformation used. 

In the present article we wish to account, even if in a model dependent way, for at 
least some microscopic mechanisms whose existence is known already. In particular we 
are interested in the rotational state wave functions for increasing angular momenta. 
It is known that an increase of the angular momentum causes an increasing effect of 
the Coriolis coupling, the latter resulting in a gradual increase in the rotation-induced 
i^-mixing. Since the presence of angular momentum introduces an extra direction in 
space, on top of the original tetrahedral symmetry, the latter is gradually more and more 
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perturbed. Studying of tetrahedral symmetry of a microscopic many-body system under 
the condition of increasing spin is a matter of a compromise between the original mean- 
field symmetry-properties and the Coriolis perturbation. As long as the Coriolis effects 
can be considered small, one may talk about the tetrahedral symmetry in the system. 

In the pure mean-field context the rotational motion has been studied in the past by 
simulating the Coriolis coupling effects with the so-called cranking term which is linear 
in angular momentum and which contains the Lagrange multiplier u in the case of the 
one- dimensional rotation. More generally, in the case of three-dimensional rotation, as in 
the present triplet of Lagrange multipliers {uxiOJy^ujz} = a; is introduced. The 

notation can be shortened to uj = un, where the unit vector n points to the direction 
of the total spin. One can show further that the term uj can be given an interpretation 
of the rotational frequency valid asymptotically for increasing spin and regular nuclear 
energy- vs. -spin dependence - where from the notation uj — )■ ajfot- 

To take into account the tetrahedral symmetry mean- field, here through /idcf, the 
pairing correlation through the term /ipair, which does not impact the symmetry consid- 
erations, and a gradual effect of the Coriolis mechanism through the {\K\ = l)-mixing 
term, h\K\=i, we finally introduce the mean-field Hamiltonian 

Kiii = ^dcf + ^pair + h\K\=U h\K\=l = "l^rot Tl - J. (9) 

As it turns out the third term in the above equation with a small c^rot? typically of the 
order of hw^ot = 0.010 MeV, is sufficient to break the time-reversal invariance and to 
introduce the (l/^-R'l = 1) i^-mixing in the wave function, which is important to obtain 



more reliably the moment of inertia 20| within the Hill- Wheeler system of equations for 
the states not far from the ground-state. It has been tested, cf. Fig. 7 in Ref. j20|, that 
the resulting spectra do not depend very much on the particular choice of values of hwrot-, 
as long as the spin values involved are not too high. Since in this paper we are interest in 
the spins of up to a dozen of the particular value of this coefficient does not play any 
essential role and we keep the above value without modification for the present purpose. 

Let us emphasise that the present use of the third term in Eq. ([9]), as compared to its 
role in the standard cranking model, is different. Whereas in the cranking-model approach 
in numerous articles on the high spin physics, Wrot plays either the role of the Lagrange 
multiplier adjusted to each new spin value and thus, on the average, increasing with spin. 



or, alternatively, is currently being used as an independent cranking variable in function 
of which observables such as e.g. single-nucleon Routhians are plotted - here it may be 
seen as a coupling constant in front of a certain phenomenological interaction term. 



B. Two-Body Model Hamiltonian for Projection Calculations 

As commented already earlier, we wish to go beyond the mean-field approximation 
to be able to take into account, at least partially, certain two-body correlations which 
have proven to be successful in a phenomenological description of not strongly-deformed 
nuclei. For this purpose, we employ the model Hamiltonian used in the Hill- Wheeler 
formalism, cf. Eqs. ([I]) — ([3]), with an auxiliary spherically-symmetric WS potential and 



separable, schematic two-body interactions as the ones employed in Ref. [20|. The two- 
body Hamiltonian of this form does not perturb the tetrahedral symmetry, whereas at 
the same time allows for including a richer structure of the nucleon-nucleon interactions. 
More specifically, we define 

H = ho + HF + Hg, (10) 

where ho is a one-body Hamiltonian composed of the kinetic energy term and the spherical 
WS potential (with the Coulomb interaction for protons). The schematic particle-hole 
(F-type) interaction. Hp, is chosen to be isoscalar and is defined by 

A>2 fj, ''"=11, p 

where : : denotes the normal ordering and r = n, p distinguishes neutrons and protons. 
Furthermore, the spatial representation of the above particle-hole type operator, F^^, is 
defined through the one-body field, 

Fj:,{r) = R^o^Y,,{9,<p), (12) 

with V^{r) and being the central part of the WS potential and its radius, respectively. 
The so-called self-consistent value [ij] of the force strength, x, common to all multipolar- 
ities, is calculated by 



X — (^n + f^p) 



where pl{r) is the density of a hypothetical spherical ground state, which is calculated 
with the filling approximation for each nucleus based on the spherical WS single-particle 
state of Hq. On the other hand, the pairing type (G-type) interaction, Hq, acts only 
within like-particles, and is given by 

Hg = -J2T.91 E^i^I,, Gil = I (14) 

T=n, p A>0 /J. ij 

where the matrix elements of the pairing type operator, G^j^, are calculated with the help 
of the standard multipole form. 



with Rn = 1.2^1/3 



The present formalism follows the main lines of Ref. |20| with a few modifications. 
Firstly, the extra one-body terms {hi = —hp — he in §3.1 in 20|) are included in Ref. 20| 
in order to cancel out the one-body exchange contributions of the multipole interactions. 
Hp + Hg- It turns out that the effect of these terms on the resultant projected spectra is 
small, so that they are neglected for simplicity. 

A slightly different deformed mean- field Hamiltonian has been used in Ref. [20[, namely 
the one derived as the Hartree approximation to the interaction (fTTll . in the form hdet = 
ho — ctA^t-^A^j, which, however, coincides with the central part of the present deformed 
WS potential only within the first order in the deformation parameters {axfj,}- We employ, 
in the present work, the general shape parameterisation based on the deformed radius in 
Eq. (jl]) with the volume-conservation condition properly taken into account. Moreover, 
the cut-off of the pairing model space is introduced directly in the operator G^J^ in Ref. 20 |: 
then one has to use the spherical single-particle energy in the cut-off function ([8]) to keep 
the spherical invariance of the Hamiltonian. We found that it sometimes causes a problem 
that the results are rather sensitive to the choice of model space for the relatively small 
pairing model space like in the present case, = Ai = 1.2 hw. In the present calculation, 
the cut-off of the pairing model space is taken into account in the step of deformed HFB 
calculation in Eqs. (jl]) — ([7]) based on the deformed WS single-particle state. Therefore, 
the cut-off function is not included explicitly in the pairing operator G^j^ anymore. 
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C. Choice of Parameters 



The deformed mean-field single-particle states are calculated in the present work using 
the Woods-Saxon potential. An often used parameterisation introduced over thirty years 



back is referred to as 'universal' (cf. Refs. 24l-l28|). We employ, in this work, a new 
improved "universal compact" set, whose parameters are listed in Table [U see Ref. 22| 
for notations. 

TABLE I: The parameters of the Woods-Saxon potential used in this work. Symbols v and vr 
refer to neutrons and protons, respectively. 





V^c [MeV] 


Kc ^'Oc [fm] 


Oc [fm] 


A 


Voso [MeV] 




J'Oso [fm] 


ttso [fm] 


V 


-52.0 


0.650 1.26 


0.64 


28.0 


49.6 





0.870 


0.70 


TT 


-53.0 


0.526 1.27 


0.71 


23.0 


49.6 





0.888 


0.86 



As for the maximum number of the harmonic oscillator shells to be used in the calcu- 
lation, we employ A^^ax = 20, which is a safe margin to accurately calculate the single- 
particle wave functions of the bound-states in the Woods-Saxon potential, and at the 

□ 

same time guarantees the convergence of the result of projection calculations [201 ]. 

The ground state deformation is determined for each nucleus by the axially symmetric 
WS-Strutinsky calculation of Ref. ^^.-^ where an algorithm allowing for the treatment 
of nuclei with weakly bound nucleons has been implemented. In the present realisation 
we calculate the strength of the seniority force, G, so as to reproduce the even-odd mass 
differences with the calculated deformation. More precisely, for the calculated equilibrium 
(ground state) deformations, and 040, we adjust the G- values in such a way that the 
calculated pairing gap, = Gr {P}) (r=n, p), agrees with the even-odd mass difference. 
Once the parameters G are fixed in this way, the usual BCS or the HFB equations (in the 
case of the cranking Hamiltonian) are solved self-consistently at any given deformation. 

As for the spherically-symmetric Hamiltonian used for the projection calculation, the 
self-consistent value x fo^' the F-type interaction in Eq. ( |T3l) is used without any modi- 
fications. We include the A = 2,3,4 components for the F-type interaction in Eq. ( fTTl) . 
because the a2o and 040 deformations are taken into account for the ground state within 
the Strutinsky method and the tetrahedral shape is described by the 032 deformation. As 
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for the G-type interaction, the monopole pairing is known to be essential. It has been 
recognised that the quadrupole pairing interaction is also important especially to describe 
the rotational motion Therefore we include the A = 0,2 components for the G-type 
interactions in Eq. ( fT^ . The strength of the monopole pairing is determined again to 
reproduce the pairing gap, i.e., = (GqI) (r=n, p), for the ground state wave func- 
tion. As for the strength of the quadrupole pairing we assume Q^^/Qq = 13.6, which is 
determined to approximately reproduce the moment of inertia of the ground state bands 
in the previous calculation 20|. Note that the values of Gr and Qq are slightly different. 



because the monopole pairing operator and the associated strength, Gr, are defined 



with respect to the deformed WS basis, while the corresponding operator GqI and its 



strength, g^, are defined with respect to the spherical WS basis. 

The model space truncation in the projection calculation is controlled by the small 
parameter e defined through the requirement, that the orbitals which in canonical rep- 
resentation have occupation probabilities vf > e are included, and, similarly, the core 
orbitals with 1 — vf < e. We have chosen e = 10~^ in the present calculation, with which 
it is confirmed that the resultant projected energies are stable within in 1 keV, which 
corresponds to six digits of accuracy for absolute energy of the present Hamiltonian in 
Eq. (fTOj) . As for the calculation of the Hill- Wheeler Eq. ([1]), the states that have smaller 
norm eigen-values than 10~^° are excluded. 



III. RESULTS OF THE CALCULATIONS 



3. fl^ 



According to the calculations in Refs. j^|9|], the tetrahedral magic numbers are Nt or Zt 
= 16, 20, 32, 40, 56-68, 70, 90-94, 112, and 136/142, and are the same for the neutrons 
and protons. Calculations which followed, Refs. j29|, figure 3, and 30[, using the universal 
Woods-Saxon mean-field Hamiltonian, suggested that, in particular, 4oZr4o and fo^rse are 
tetrahedrally-symmetric in their ground-states, whereas tetrahedral minima lie about 1 
MeV above the ground-state in ^^J^Ztiq. Skyrme Hartree-Fock calculations for the latter 
nucleus predict the possibility of the tetrahedral minima being the lowest - depending 



strongly on the choice of the model parameters 13|. More recent 



the calculations of the tetrahedral spectra in loS'ii^Zr, in Ref. |3l|- The predicted shape 



y, we have reported on 
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coexistence in the Zirconium region which includes tetrahedral symmetry minima may 
give rise to, among others, the presence of isomeric states. Quite recently the experiment 
have been performed for unstable nucleus ^'^^Zr and the results are compatible with the 



existence of an isomeric state 32[ whose nature is being debated. 



A. Remarks about Symmetry Properties and Quantum Rotors 

One of our goals is to determine whether the microscopic calculations which combine 
various advanced techniques of the nuclear quantum mechanics reproduce the excita- 
tion pattern predicted by group-representation theory. Because all of the two-fold and 
four-fold degenerate single-particle states are occupied with equal probabilities for the 
doubly-closed shell tetrahedral configurations, the totally symmetric, so-called tensor Ai 
irreducible representation of the tetrahedral point-group (as opposed to the spinor 
irreducible representations characterising the symmetry properties of the single-nucleonic 
wave-functions within the tetrahedral double-point group T^, cf. also Table VIII in Ap- 
pendix) can be expected as the resulting symmetry of the full system in its lowest rota- 
tional band. The situation remains the same if the seniority-type pairing interaction is 
effective for fully paired even-even nuclei. 

An ideal tetrahedral (T^ -symmetric) classical rotor is often referred to as 'spherical', 
because its moment of inertia tensor is diagonal with all elements strictly equal. In the 
case of a quantum rotor, the notion of the inertia tensor cannot be strictly-speaking 
defined since the only quantum observables directly associated with rotational motion of 
such an object are the energy and angular momentum. The specific spectral properties of 
quantum rotors with point- group symmetries have been actively studied in relation to the 
TetraNuc Collaboration activities in recent years. For instance, examples of the octupole- 



symmetric quantum-rotor spectra have been presented in Ref. 33|; simf 



specifically tetrahedral-symmetric quantum rotors can be found in Ref. 34] whereas the 



ar examples for 



underlying tensor formalism and the general form of the reduced matrix elements are 



presented in Ref. 35|. Observe that in contrast to the 'usual' quantum rotor Hamiltonians 
discussed in the literature, which are quadratic forms of the angular momentum operators, 
the tetrahedral (or other octupole-symmetric) rotor Hamiltonians are specific third-order 



12 



forms expressed in terms of the operators {Ix,Iy,Iz} (equivalently of {J_i, Jq, J+i} using 
spherical tensor representation). Furthermore, relations between the energy spectra of 
the quantum rotors and the associated properties of classical rotors have been discussed 



in Ref . 



36| 



It can be shown that the spectra of tetrahedral-symmetric structure-less quantum rotor 
are composed of the (2/ + l)-degenerate states for each given spin /, cf. Ref. 3J]. In this 
sense the spherical-symmetry of the classical tetrahedral rotor mentioned above and the 
symmetry of the tetrahedral quantum rotor can be seen as analogous. According to group- 
theory considerations, each of the (2/+l)-degenerate states of the -symmetric structure- 
less rotor of any given / belongs to a certain specific irreducible representation of the group 
in question. Among five of those irreducible representations, the Ai (scalar) representation 
contains states with the following characteristic set of spin-parity combinations [cf. Table 
VI in the Appendix (also e.g. Ref. 
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39|)]: 



0+, 3", 4+, 6+, 6", 7", 8+, 9+, 9", 10+, 10", 11", 2 x 12+, 12" 



(16) 



As it is discussed in more detail below, such characteristic spectra are indeed realised for 



3l|. 



the lowest energy band in the results of our microscopic calculation, cf. also Ref. 

It is worthwhile mentioning that there exist certain extra discrete symmetries for the 
tetrahedral shape nuclear mean-field configurations applying within the cranking model. 
They are referred to as doublex and triplex (as opposed to the 'usual' simplex- symmetry 
applying to the cranking model for the pear-shape symmetric nuclei). The corresponding 
quantum numbers are useful to further classify the characteristic spectra as discussed in 
Refs. 
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4l|, but this issue goes beyond the scope of the present article. 



In this work we report on the more extensive and detailed investigations for three 
doubly-closed tetrahedral-shell configurations in nuclei: ^^°Yb {Z = 70 and = 90), 
the already mentioned ^^^Zi, as well as the heavier tetrahedral-symmetric nucleus in the 
Actinide region, ^^^Th (Z = 90 and = 136). The axially symmetric octupole (aso) 
deformed states of the latter nucleus have also been studied in Ref. 20|]. 
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B. Tetrahedral-Symmetry in °°Yb: Phonon-vs.-Rotation-Like Structures 



Following the procedure of Sec. Ill Ct we obtain the following parameters for the pair- 
ing force and residual interaction given in Table HIl where the calculated deformation 
parameters for the ground state are also shown. 

020 "40 [MeV] Z\p [MeV] [MeV] Gp [MeV] x [MeV-^] [MeV] gl [MeV] 
0.194 0.031 1.265 1.370 0.1512 0.1732 2.643x10"^ 0.1468 0.1675 

TABLE II: The calculated ground state deformation parameters (a20) ct4o)) the 4-th order even- 
odd mass difference (/^nj ^p); and the force strength parameters determined based on them for 
^ygYbgg. Thc valuc ~ -'-^•6 taken for the ratio of the quadrupole and monopole pairing. 




Let us begin by presenting the results of our projection calculations for the ground 
state band in ^^°Yb; the corresponding results are shown in Fig. [H where comparison 
with the experimental data can also be found. In this calculation, we use the HFB type 
wave function with a small Coriolis (|AJr| = l)-coupling /iWrot = 0.010 MeV parameter, 
as discussed in Sec. Ill A[ We choose the cranking axis as the y-axis, which is perpendicular 
to the symmetry axis (^-axis). The agreement of the calculated spectra with the experi- 
mental data is acceptable but discrepancies as compared to the experiment increase with 
spin. In particular, the observed moment of inertia increases as a function of spin, but the 
calculated moment of inertia is fairly constant. This trend was already found for other 
nuclei in Ref. 20| and may partly reflect the fact that in contrast to the cranking model, 
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within which the pairing correlations systematically decrease with cranking frequency and 
thus increase the moment of inertia, in the present model with projection from one mean- 
field state the mentioned mechanism does not exist. In the present calculation of the 
ground state the axial symmetry is broken only by a small (|AJi'| = l)-mixing term, and 
the mixing effect of the K quantum number in the Hill- Wheeler Eq. (jj]) remains small. 
Since this particular aspect is of secondary importance of the present project we accepted 
the disagreement in question as remaining under control but without consequences for the 
main conclusions. Doing so we may profit from a technical advantage: It is sufficient to 
use relatively small numbers of points for Gauss quadratures with respect to the Euler an- 
gles when performing the angular-momentum projection calculation, e.g., = = 16 
and Nfj = 50. 

We proceed to examining the tetrahedral nuclear configuration and related rotational 
states. The simplest way to construct the tetrahedral shape in the surface parameterisa- 
tion in Eq. (jl]) is to set 032 = — 0^3-2 as the only non-zero deformation parameters. 
In the coordinate system chosen the upper [z > 0) and lower [z < 0) sides of the tetrahe- 
dron are parallel to the y- and respectively, and the 2-axis is along the line that 
connects the middle points of these two facing sides. Again, we have chosen the y-axis 
for cranking with a small frequency hui^ot = 0.01 MeV. Obviously, in this case the axial 
symmetry is strongly broken depending on 032 and thus higher order quadratures in the 
projection calculations are necessary. We take = N-y = Np = 64 for ^32 = 0.1 — 0.2, 
Na = N^ = Ni3 = 84 for ^32 = 0.25-0.30, and N„ = = N/i = 104 for 032 = 0.35-0.40 
for the calculation in this nucleus. 

In Figure |2] we illustrate the result of excitation energies of the projected eigen-states, 
cf. Eq. ([1]), whose norm is not very small, i.e. not less than 10~^ of that of the ground 
state. They compose the lowest-energy sequence, and calculated for the pure tetrahedrally 
deformed configuration with deformations 032 = 0.05, 0.10 and 0.20. As it is seen from the 
figure, only the specific spin-parity combinations appear, cf. Eq. (fT6l) . for all deformations. 
This pattern is characteristic for the tetrahedral symmetry and is quite different from the 
one of the usual quadrupole deformation. Indeed, the spectrum is composed of states 
characterised by O"*", 3~, 4"'', 6^, the states with J = 1, 2 and 5 are missing in the 
figure since their norms are too small and/or they lie much higher in energy as compared 
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FIG. 2: Examples of calculated spectra of tetrahedral states belonging to the lowest energy part 
in the spectrum in ^^"Yb with 032 = 0.05 (left), 032 = 0.10 (middle) and a32 = 0.20 (right), test- 
ing the dependence of the excitation energies on the tetrahedral deformation [cf. also Eq. ()16p ]. 
The positive (negative) parity states are denoted by the solid (dotted) lines. Note the difference 
of ordinate-scales between the left and right panels. 
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to the discussed lowest-energy sequence. 

A very small tetrahedral deformation of ci;32 = 0.05 corresponds to nearly spherical 
form. The existence of collective excitations in spherically symmetric nuclei has long been 
recognised in terms of the vibrational modes, which results in an equidistant spectrum 
composed of multiplet of states, within the simplest harmonic representation, in terms of 
the vibration-quanta: The phonons, see e.g. Sect. (6.3.2) in Ref. (37|. 

The spectrum with smallest deformation 032 = 0.05 in Fig. [2] is more vibrational-like, 
while that with 032 = 0.20 is approaching to the rotational-like spectrum; the one with 
^32 = 0.10 is in-between. In fact, (4+,6+), (6", 7^,9,-), (8+, 9+, 10+, 12+) ... states in 
Fig. |2] a) having the same parity can be grouped together, and would be interpreted 
as slightly perturbed two-phonon, three-phonon, four-phonon ... multiplet structures, 
respectively, of an elementary mode of the 3~ vibrational excitation. Moreover, the exci- 
tation energies of 3~, 6+, 9~, 12+ form a rather linear dependence as a function of spin; 
the dependence resembles the pattern expected for the multi-phonon excitation. Note, 
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however, that only the specific spin states among the multi-phonon muhiplets appear, 
which is a consequence of the tetrahedral symmetry. In contrast in Fig. |2]c), the states 
with the same spin value are nearly degenerate, which is a specific feature of the ideal 
rotor, and, at the same time, approximately follow the quadratic energy-vs.-spin relation, 
oc/(/ + l). 

Let us emphasise, that both the parity and the angular-momentum projections were 
essential for obtaining the tetrahedral-symmetry pattern predicted by the group theory. 
This symmetry pattern seems rather typical for the present model Hamiltonian: We obtain 



31| and 226r]n]^ below. In the present 



similar pattern also for other nuclei, e.g., ^'^^'^^^Zr 
work, we concentrate on the lowest energy sequence and we do not enter the discussion 
of the group theory aspects. Instead let us only mention that representations other than 
Ai must be expected for excited bands; also - one may expect, that the symmetry in the 
case of the non-doubly-closed shell nuclei could be manifested less strongly. 
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FIG. 3: Comparison of the spectra of tetrahedral states with (right) and without (left) the 
number projection (NP), which are shown in the same way as in Fig. [2] but with Q32 = 0.15. 



In Figure [3] we compare the results of tetrahedral spectra with and without the particle- 
number projection (NP) related to pairing formalism. Although the moments of inertia 
(the slopes) are slightly different, the characteristic properties of the spectra are exactly 
the same in the two calculations. We conclude that the effect of the particle number 
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projection is small and we do not apply it in the rest of the article. 



Transition to Ideal Rotor and Moments of Inertia in ^^°Yb 
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FIG. 4: Calculated spectra of tetrahedral states in ^^'^Yh with 032 = 0.10, 0.15, 0.20, 0.25, 0.30, 
and 0.35, respectively, for a), b), c), d), e), and f). The dotted line in each panel denotes an ideal 
/(/+1) sequence going through the first excited 3~ state. Note that almost exact degeneracies for 
/ = (6+,6~), (9+,9-), (10+, 10"), (2x12+, 12~) states are obtained for 032 > 0.25 demonstrating 
the nearly perfect rotor character of the rotational excitation of the system. 



The fact that the tensor of inertia of an ideal classical tetrahedral rotor is diagonal with 
all components equal ('spherical rotor') suggests that the tilting direction of the cranking 
axis in Eq. iQ may not affect the tetrahedral spectra, at least to the extent in which the 
Coriolis alignment effects can be neglected, i.e., for not too high spins. In order to test 
this conjuncture, we investigated the projected spectra from the cranked mean-field state 
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with the different tilted cranlcing axes. We have varied the cranlcing axis in our coordinate 
system, i.e. the vector n in the term, h\K\=i = — i^rot n ■ J, in Eq. is changed by 



n = (sin6'sinv9,cos6',sin6'cosv3), < < 90°, < v? < 45°. 



:i7) 



We found that the differences of the resulting spectra for the lowest energy sequence are 
negligible within the accuracy of our calculation; i.e. the nature of 'spherical rotor' is nu- 
merically confirmed. More generally, the projected spectra for the tetrahedral symmetric 
nuclei do not depend on the \AK\ = 1 coupling term, both the strength u^t and and the 
direction n of the tilted cranking long as cjrot is small. 
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FIG. 5: Calculated spectra of quadrupole deformed states in ^^°Yb with 020 = 0.10, 0.15, 0.20, 
0.25, 0.30, and 0.35, respectively, for a), b), c), d), e), and f). The dotted lines in each panel 
denote the ideal I {I + 1) spectra going through the first excited 2"'" state. 



In Figure H] we show the calculated excitation energies for selected values of the tetra- 
hedral deformation. The ideal rotor spectra with the energies proportional to /(/ -|- 1) 
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and containing the calculated first excited 3~ states are also shown by the dotted lines. 
This figure clearly shows that the spectra exhibit the gradual transition from linear to 
parabolic spin dependence with increasing the tetrahedral deformation: The almost ideal 
rotor spectrum is realised for > 0.25. 

Now we compare the tetrahedrally-symmetric spectra with those of the quadrupole 
deformation. In Figure |5] we show the results of calculated spectra obtained by the 
angular-momentum projection from the pure quadruple deformed states, where all the 
deformation parameters are set to zero except 0:20 (no parity projection is required in this 
case). The projection calculation tends to give good rotational spectra, but the result 
with small deformation, 020 = 0.10, considerably deviates from the one for the pure 
rotor spectra. Thus the gradual transition from the linear to parabolic energy-vs.-spin 
dependence is seen also for the calculation of the quadrupole deformation. 

Although the moment of inertia is not any quantum-mechanical observable, certain 
quasi-classical analogies often found in the literature allow to define and estimate the 
corresponding values. Here we define this parameter, J', through 

m-EiO) = ^-^^. (18) 

It is well-known that the moments of inertia of observed rotational band near the ground 
state are about (or even smaller than) half of the classical rigid-body value. This large 

n n 

reduction is supposed to be due to the pairing correlations [1|, |42[. In fact, the moments 
of inertia extracted from the high-spin states, where the pairing correlations are believed 
to be quenched, are known to be close to the rigid-body value, although some deviations 



attributed shell effects exists, see e.g. Ref. 43|. Therefore, it is instructive to investigate 



the moment of inertia in the case of the tetrahedral rotor. In Figure the moments of 
inertia calculated from Eq. ( IT8l) are plotted as functions of the tetrahedral (left) and the 
quadrupole (right) deformations, where they are estimated from the calculated 3~ and 2+ 
excitation energies, respectively. The results with neglecting the pairing correlation and 
the rigid-body value are also included. An irregular behaviour for the unpaired (A = 0) 
moments of inertia, i.e., at 020 ~ 0.20 — 0.25, is due to the fact that the level crossings 
near the Fermi surface occur. In order to illustrate the possible correlation between the 
moments of inertia and the intensity of pairing correlations measured with the help of 
the pairing gaps, the calculated pairing gaps at corresponding deformations are shown in 
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FIG. 6: Moment of inertia estimated from the calculated spectra for the pure tetrahedral 
states (left) and the pure quadrupole states (right) as functions of the deformation parameters 
in ^^"^Yb. The energy of the first excited 3~ (2"^) is used for estimation of the former (latter). 
The results with the pairing correlations artificially set to zero are also included. The classical 
rigid-body moments of inertia in function of the deformation parameters are shown as solid lines. 
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FIG. 7: The self-consistent neutron and proton pairing gaps for the tetrahedral states (left) and 
the quadrupole states (right) as functions of the respective deformation parameters in ^^^Yb. 



As it becomes clear from Fig. [6l the moment of inertia increases rapidly with increasing 
deformation for both the tetrahedral and quadrupole shapes, indicating that the picture 
of the good rotor emerges for larger deformation. Observe that the values of the moments 
of inertia at the two considered shape configurations are rather similar when the pairing 
correlations are included, and are both much smaller than the rigid-body values even at 
the largest value of the deformation parameters. If the pairing correlations are set to 
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zero, the calculated moment of inertia for the quadrupole shape becomes much larger and 
approaches to the rigid-body value as it is observed from the high-spin limit. However, 
the effect of pairing correlation on the inertia for the tetrahedral shape is rather small, 
even though the pairing gap takes more or less the same values as in the case of the 
quadrupole shape. The reason why the moment of inertia is small and is affected very 
weakly by the pairing correlation may be because the chosen nucleus in this case is the 
tetrahedral doubly-closed shell nucleus. The shell gap is ~ 1.5 — 2 MeV and is larger than 
the pairing gap, which is in contrast to the case of quadrupole deformation, where the 
mean single-particle level spacing is much smaller than the pairing gap. 



D. Tetrahedral spectra in ^Zr 



The nuclear potential energy surfaces for the doubly-magic Zirconium nuclei have been 



studied in Ref. 



and the corresponding illustrations obtained using the phenomenolog- 



ical approach with the Woods-Saxon mean-field Hamiltonian can be found in figure 3 of 



the above reference, 
can be found in Ref. 



he symmetry-oriented discussion of the corresponding shell-effects 



44l |. A discussion of the static-energy properties in a few nuclei in 



the vicinity of ^^°Zr using Hartree-Fock approach can be found in Ref. 45|, whereas the 
tetrahedral rotational properties, specifically for the nucleus ^^°Zr, have been studied us- 
ing the cranking-Skyrme-Hartree-Fock method in Ref. 46|] and using the methods similar 
to that of the present article in Ref. 31|. 



In the present work the method of calculation is essentially the same as in [3l|, except 
that the different Woods-Saxon Hamiltonian parameter set is used. The mean-field pa- 
rameters and the force strengths used in the present calculation are given in Table UTTl The 
calculated pairing gaps are used because the experimental even-odd mass differences are 
not available for this unstable nucleus. The numbers of nodes for the Gaussian quadra- 
tures are chosen to be = = = 64 after verifying the stability conditions for the 
final result. 

In Figure [8] we show the excitation energies for various deformations in function of 
angular momentum, obtained using the angular-momentum and parity projection tech- 
niques. The cranking axis is chosen to be the y-axis [6 = 0° in Eqs. and ( |T7I) ]. As it is 
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020 "40 [MeV] Ap [MeV] Gn [MeV] Gp [MeV] x [MeV"^] ffg [MeV] c/{; [MeV] 



0.333 -0.026 1.129 1.113 0.1471 0.2698 4.785x10"^ 0.1407 0.2625 



TABLE III: The calculated ground-state mean-field parameters {020, 040, ^p)) ^-nd the 
force strength parameters determined based on them for ^IqZtjq. The value S'2/5'0 ~ ^"^-^ 
taken for the ratio of the quadrupole and monopole pairing. 
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FIG. 8: Calculated spectra of tetrahedral states in ^^°Zr with 032 = 0.10, 0.15, 0.20, 0.25, 0.30, 
and 0.35, respectively, for a), b), c), d), e), and f). The figure is similar to that in Ref. [31 1 



but only the lowest band is selected and the results for larger deformation are included. Note 
that almost exact degeneracies for / = (6"'", 6~), (9"'", 9~), (10+, 10^), (2 x 12+, 12^) states are 
obtained for 032 > 0.30. 



seen from the Figure, the over-all pattern of the excitation scheme resembles the one in 
16OYI3 However, compared with the results for ^^°Yb of Fig. HI the transition to the ideal 
rotor occurs slower, i.e., it occurs at the larger deformation in the lighter system ^^°Zr. 
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FIG. 9: Moment of inertia estimated from the calculated spectra for the pure tetrahedral states 
(left) and the pure quadrupole states (right) in ^^^Zi. The energy of the first excited 3~ (2"*^) is 
used for estimation of the former (latter). The rigid-body moment of inertia are shown as solid 
lines. 



In particular, the energy-vs.-spin relation resembles a rigid rotation only for > 0.30. 
Compared with the results in Ref. 3l|, those in the present work are very similar, in- 
dicating that the choice of the parameter set of the WS potentials very little affects the 
rotational properties of a tetrahedral symmetric nucleus - provided a realistic choice of 
parameters is used. We have checked the dependence of the projected energies on the 
tilting angle of the cranking axis also for the nucleus ^^°Zr (cf. also figure 1 in Ref. j4l|). 
Again, the result is found to stay the same, when the tilting angle is changed in the same 
way as in i^°Yb [cf. Eq. ffTTl) ]. 

In Fig. ini the moments of inertia for the tetrahedral and quadrupole deformations are 
compared for the case of ^^^Zt; only the results with pairing correlations are presented. As 
it is seen, the values of the moments of inertia for two types of deformations are slightly 
different; the ratio i7(3~)/i7rigid for the tetrahedral shape is considerably smaller than 
the ratio Jl{2'^)/ J^rigid for the quadrupole shape. One of the reasons may be traced back 
to the somewhat larger shell gap at Z = 40 for the tetrahedral shape, so that the ratio 
j7'(3~)/i7rigid in ^^"^Zr is smaller than that in ^^°Yb. Furthermore, the pairing gaps for the 
quadrupole shape with 020 > 0.20 are somewhat reduced in ^^°Zr, which increases the 
quadrupole moment of inertia; those combined effects for the two types of shapes may 
make the difference of their behaviour in ^^°Zr compared to the case of ^^°Yb (and of 
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Th, see below). 



E. Tetrahedral spectra in ^'"Th 

The calculation procedure for ^^^Th is the same as that for ^^^Yh and ^^°Zr. The 
mean-field single-nucleon energies for this particular nucleus in function of tetrahedral 
deformation can be found in figure 4 of Ref. {4^. The parameters determined by the 
Woods-Saxon-Strutinsky calculation, and those of the force strengths are tabulated in 
Table ITVl They are similar to the ones used in Ref. {20], where the a^o deformation was 
also taken into account. Slightly different values of the parameters as compared to those in 
Ref. 20j are mainly due to the fact that the different parameter set of the WS potential is 
employed. The computing time of the projection calculation increases dramatically with 
increase of the nucleon number as well as the numbers of the mesh points for Gaussian 
quadratures. We have carefully tuned the latter numbers for '^'^^Th to obtain the same 
accuracy as in the case of ^^^Yh. Thus, we take = = 64 and Nj^ = 74 for 
^32 = 0.10 - 0.15, Na = = 104 and iV^ = 84 for 032 = 0.20, = = 124 and 
iV^ = 104 for = 0.25 - 0.30, and = = = 124 for 032 = 0.35 - 0.40. 

020 [MeV] Ap [MeV] Gn [MeV] Gp [MeV] x [MeV"^] [MeV] [MeV] 

0.161 0.093 0.814 0.830 0.09772 0.1289 1.744x10""^ 0.09568 0.1267 



TABLE IV: The calculated ground state deformation parameters {020, 040)5 the 4-th order even- 
odd mass difference (Z\n, ^p); and the force strength parameters determined based on them for 
^ggTh^gg. Thc valuc g^/ 9q = 13.6 is taken for the ratio of the quadrupole and monopole pairing. 

The result of our angular-momentum and parity projection calculations for various 
tetrahedral deformations are shown in Fig. [TOl The characteristic features of the spectra 
resemble those of ^^^Yb and ^^'^Zr. Again, the energy-vs.-spin dependence has approxi- 
mately linear behaviour for smaller deformations, whereas it approaches a parabolic form 
at increasing deformation. Comparing the energy-vs.-spin dependence of these three nu- 
clei, the transition from a linear to parabolic spin dependence occurs at smaller defor- 
mation in nuclei with larger mass number. More precisely, the energy-vs.-spin depen- 
dence becomes almost parabolic in the following proportions: at 032 ~ 0.15 in ^^^Th, at 

25 



> 



6f 
I 



6 
4 





> 6 







a) 


□ 


"In 


• 

• 








••■ • 032 = 0.10 




n=+ • 




;:=- □ 




, 1(1+1) , 



10 12 







c) 








032 


= 0.20 


..1 


.a' 







6 

I[h] 



10 12 




d) 






032 = 0.25 




..a 




«■-■■■ 









6 



10 



12 



FIG. 10: Calculated spectra of tetrahedral-symmetry states in '^'^^Th with 032 = 0.10, 0.15, 
0.20, and 0.25 respectively, for a), b), c), and d). Note that almost exact degeneracies for 
/ = (6+,6-),(9+,9~),(10+,10-),(2 X 12+, 12") states are obtained for 032 > 0.20. 



032 ~ 0.20 in ^^'^Yb, and at 032 ~ 0.25 in ^^^Zt. This is intuitively acceptable because 
the concept of the symmetry-breaking is more and more appropriate for heavier nuclear 
systems. 
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FIG. 11: Moments of inertia estimated from the calculated spectra for the pure tetrahedral 
states (left) and the pure quadrupole states (right) in ^^^Th. The energy of the first excited 3~ 
(2+) is used for estimation of the former (latter). The classical rigid-body moments of inertia 
are shown as solid lines. 
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The calculated moments of inertia as functions of quadrupole and tetrahedral defor- 
mation parameter are illustrated in Fig. [TTl where only the result including the pairing 
correlation is shown. Again, the moments of inertia with the pairing correlation are con- 
siderably smaller than the rigid-body values and they increases with deformation. The 
values of moment of inertia are rather similar for the tetrahedral and quadrupole shapes 
in ^^^Th as in the case of ^^°Yb in Fig. El 



IV. SUMMARY 



We have studied the rotational nuclear properties for the pure tetrahedral deformation 
by the angular-momentum and parity projection method employing the realistic Woods- 
Saxon mean-field potential and the schematic separable two-body interaction consistent 
with it. In this work we have chosen the tetrahedral doubly-closed shell nuclei ^^°Yb 
{Z = 70 and = 90), ^^^Zr {Z = 40 and A^ = 70) (see also Ref. [aj for this nucleus) 
and ^^^Th {Z = 90 and A^ = 136) as illustrative examples. We have found out that the 
characteristic spectra for the totally symmetric representation, i.e. the Ai irreducible rep- 
resentation of the tetrahedral group T^, are obtained with the specific sequence composed 
of O"*", 3~, 4"*", 6~^, 6~, 7", 8"^,- ■ ■ . These spectra are consistent with those of the simple 
tetrahedral rotor, although the spectra are more vibrational- like for small deformations. 
However, it is important to emphasise that only the specific spin-parity combinations 
appear in the projection calculations in this work in agreement with the group-theory 
predictions, approaching the rotational pattern closer and closer with increasing tetrahe- 
dral deformation. 

The Coriolis |AA'| = 1 mixing is introduced in order to break the time-reversal in- 
variance and to obtain more reliable estimate of moment of inertia; only the slope of 
energy-vs.-spin relation changes and the qualitative features are not affected by this mix- 
ing. It has been also checked that the results of projected spectra are independent of the 
tilting angle of the cranking axis, which is consistent with the picture that the tetrahedral 
rotor possesses a 'spherically-symmetric tensor of inertia'. 

Stable rigid-rotor energy dependence [oc J(J + 1)] appears for larger deformations in 
the range studied in this article, but an approximately linear dependence is obtained with 



27 



decreasing deformations. This transition between the vibrational- hke and the rotational 
pattern occurs at smaller deformation for nuclei with larger mass number. The moment 
of inertia for the tetrahedral shape increases with increasing deformation, but it is much 
smaller than the rigid-body value irrespective of the pairing correlations. The impact of 
the pairing seems rather limited for tetrahedral deformed nuclei at least for the doubly- 
closed shell nuclei, with relatively large single-particle energy-gaps. 

In this work, we have concentrated on the lowest-energy rotational sequences and 
shown by the calculations that their properties resemble the ones characteristic for the 
totally symmetric, the so-called Ai irreducible representation of the tetrahedral point- 
group, Td, believed to be characteristic for the lowest-energy sequence of states in fully 
paired even-even nuclei. From the group theoretical consideration, it is expected that the 
other irreducible representations, i.e., A2, E and/or Fi, F2, would appear in the excited 



311 ]. or in the spectra for 



rotational bands of even-even nuclei, see the results of Ref. 
odd-odd nuclei. Moreover, there exist different types of representations associated with 
the so-called double tetrahedral group, specific for the odd nuclei, see Appendix. 

We believe that both the experimental analysis as well as the theoretical calculations 
of the discussed spectral properties will need to include more levels in the future. For 
this purpose the remaining irreducible representations of the symmetry group may need 
to be studied compared to the scalar representations that we focused on in this article. 
The analysis of the irreducible representation structure of the solutions is important for 
the next step of the analysis, which would consist in calculating the electromagnetic 
transitions and the branching ratios, the observables which change rapidly with symmetry 
of the system. This step of the analysis will be essential for establishing the experimental 
criteria of determining the presence of tetrahedral symmetry in the physics of subatomic 
systems. 
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Appendix: Spin-parity relations in a tetrahedrally symmetric rotor 



Although it may be considered a textbook matter, some spin-parity properties of the 
rotational energies of the tetrahedrally symmetric rotor will be summarised in this Ap- 
pendix, to facilitate the comparison between the results of the microscopic calculations 
with the projection techniques as obtained in this article and the group-theory expecta- 



tions (see e.g. Ref. 48|). The representation of the rotor states with definite spin-parity 
Jvr (tt = ±), D^^^\ which have a certain symmetry governed by a group G, can be de- 
composed into its irreducible representations, (z = 1, ■ ■ ■ , M), with multiplicity a-^'^^; 

M 



(A.l) 



i=l 



The multiplicity can be calculated by the standard formula [48|, 



(/tt) 



1 1 ^ 



(A.2) 

where the number A^^ is the order of the group G, Xin{R) and Xi{R) a-^e the characters of 
the representations D^^'^^ and Di, respectively, for the group element R, and the quantity 
Qa denotes the number of elements in the class a, whose representative element is Ra- Note 
that the decomposition (lA.ll) is performed by a unitary transformation in the (2/ + 1) 
dimensional space of the rotor wave functions for given Jvr; more precisely, a specific 
combination of the A'-mixing generates each irreducible representation. 



Td 


E 


cm 


^2(3) 




54(6) 




1 


1 


1 


1 


1 


A2 


1 


1 


1 


-1 


-1 


E 


2 


-1 


2 








Fi{Ti) 


3 





-1 


-1 


1 


F2{T2) 


3 





-1 


1 


-1 



TABLE V: Character table for the tetrahedral group T^- Taken from Ref. 
C2 = si and Fi^2 are sometimes denoted as Ti_2)- 



48|] (note that 
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/+ 


n+ 


1+ 


2+ 


3+ 4+ 5+ 


6+ 


7+ 


8+ 




in+ 


11+ 


12+ 13+ 14+ 15+ 16+ 


^1 


I 


n 


n 





1 





\ 


n 


1 




1 


n 


2 


1 


1 


1 2 


A2 











1 








1 


1 





1 


1 


1 


1 


1 


1 


2 1 


E 








1 





1 


1 


1 


1 


2 


1 


2 


2 


2 


2 


3 


2 3 







1 





1 


1 


2 


1 


2 


2 


3 


2 


3 


3 


4 


3 


4 4 


F2{T2) 








1 


1 


1 


1 


2 


2 


2 


2 


3 


3 


3 


3 


4 


4 4 


/- 




\- 


2" 


3- 


4- 


5- 


a- 


7- 




tj 




11- 


12" 


13- 


14- 


15- 16" 


A^ 
^1 











1 








1 


1 





1 


1 


1 


1 


1 


1 


2 1 


Ao 


1 











1 





1 





1 


1 


1 





2 


1 


1 


1 2 


E 








1 





1 


1 


1 


1 


2 


1 


2 


2 


2 


2 


3 


2 3 


Fi{Ti) 








1 


1 


1 


1 


2 


2 


2 


2 


3 


3 


3 


3 


4 


4 4 


F2{T2) 





1 





1 


1 


2 


1 


2 


2 


3 


2 


3 


3 


4 


3 


4 4 



TABLE VI: The number of states af'^^ belonging to the five irreducible representations of 
for integer spins; those for each parity are separately shown. 



The group has five irreducible representations and classes, whose representative 
elements are E, C2 (= 5*1), C3, a^, and 5*4; see Ref. 



48 



for the notation. The characters 



for each irreducible representation are listed in Table |V] for completeness, and those for 
the rotor representation 49| are as follows; 



sm 



(2J+l)7r 



K= 



sm 



(A.3) 



X/7r(0-d) = TT X X/7r(C2), X1-k{S ^ = T\ ^ XlA^ ^ ■ (A.4) 

Combination of the characters in Eqs. flA.3|) — flA.4p and in Table IVl with the formula flA.2|) 
leads to the multiplicities, ai"^^ ^ which are summarised in Table IVl] for integer spins up 



to / = 16. It is easy to verify that a)^^" 



a 



ci~p^\ In the 



table, a^^^"* = means that such states are not allowed, and a^"' = 2 means that states 
are doubly degenerate. In this way the characteristic spin-parity for the Ai representation 
in Eq. f|T6|) follows. 
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a 




E 


(^3(8) 


^2(3) 




54(6) 




2 


-2 


1 -1 








^/2 -V2 




2 


-2 


1 -1 








-^2 \/2 




4 


-4 


-1 1 












TABLE VII: Character table specific for the extended (also called 'double') tetrahedral group 
TJ^. The second entry in the corresponding columns denotes the characters of extended elements. 



[Taken from Ref. 



39( 1 (note that -©1/2) -E'5/2 ^"^^ ^^3/2 ^-I's sometimes denoted as E'^, E'2 and G') 



1+ 


1 + 
2 


3 + 
2 


5 + 
2 


7 + 
2 


9 + 
2 


11 + 
2 


13 + 
2 


15 + 
2 


17 + 
2 


19 + 
2 


21 + 
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23 + 
2 


25 + 
2 


27 + 
2 


29 + 
2 


31 + 
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Ei/2{E'i) 


1 
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1 


1 


1 


1 


2 


2 


1 


2 


3 


2 


2 


3 


^5/21 E'2) 
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1 
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2 
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1 
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2 


2 


2 


2 
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3 
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2 


2 
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3 


3 
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5 
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1- 
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7- 
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9- 
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11- 
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13- 
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2 


21- 
2 


23- 
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2 


27- 
2 


29- 
2 


31- 
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1 
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2 
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2 


2 
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2 


2 


3 


3 


Eh/2{E2) 


1 








1 


1 


1 


1 


1 
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2 
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2 


3 


2 


2 


3 


Gz/2{G') 





1 


1 


1 


2 


2 


2 


3 


3 


3 


4 


4 


4 


5 


5 


5 



TABLE VIII: The number of states a^^'^^ belonging to the three irreducible representations 
specific for T? for half-odd integer spins; those for each parity are separately shown. 



As for the excitations of odd nuclei, i.e., for half-odd integer spins /, the same cal- 
culation can be done, but one has to consider the extended rotation group (th e double 
group) , see e.g. Ref. [39| (or equivalently, the two- valued representations |48|). In the 
extended group, the number of elements is doubled by extending the range of rotational 
angle about an axis from 27r to 47r, because the 2tt rotation is not the identity operation 
but changes sign for the rotor states with half-integer spins. The character table specific 
for the double tetrahedral group Tj^ is shown in Table IVlTl and the resultant multiplicities 
are given in Table IVIIIt where a^^"* = a^^'' and o'q'^^ = can be easily confirmed. 

1/2 5/2 3/2 3/2 

Table IVIIII can be also used to see how a spherical single-particle orbit j'^ decomposes 



31 



into the two-fold (-^1/2, 5/2) and four-fold (6*3/2) degenerate orbits for finite tetrahedral 
deformation. 
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